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Abstract 

The measurements of pulsar frequency second derivatives have shown that they are 10 2 — 10 6 times larger than 
expected for standard pulsar spin-down law. Moreover, the second derivatives as well as braking indices are even 
negative for about half the pulsars. We explain these paradoxical results on the basis of the statistical analysis of 
the rotational parameters v, v and v of the subset of 295 pulsars taken mostly from the ATNF database. We have 
found a strong correlation between v and v for both v > (correlation coefficient r ~ 0.9) and v < (r ~ 0.85), as 
well as between v and v (r « 0.6 -j- 0.7). We interpret these dependencies as evolutionary ones due to v being nearly 
proportional to the pulsars' age. 

The derived statistical relations as well as "anomalous" values of v are well described by assuming the existence 
of long-time variations of the spin-down rate. The pulsar frequency evolution, therefore, consists of secular change of 
v ev (t), Vev(t) and v £V {i) according to the power law with n ~ 5, the irregularities, observed within the timespan as 
timing noise, and the non-monotonous variations on the timescale of several tens of years, which is larger than that 
of the timespan. It is possible that the nature of long-term variations is similar to that of short-term ones. The idea 
of non-constant secular pulsars' braking index n is also analysed. 
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1. Introduction 

The spin-down of radio pulsars is caused by the 
conversion of their rotation energy into emission. 
According to the "classical" approach, their rota- 
tional frequencies v evolve obeying the spin-down 
law v — —Kv n , where K is a positive constant that 
depends on the magnetic dipole moment and the 
moment of inertia of the neutron star, and n is the 
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braking index. The latter can be determined ob- 
servationally from measurements of v and v as 
n = vvfi i . For a simple vacuum dipole model of 



the pulsar magnetosphere n = 
decreases this value to n = 1; 



3; the pulsar wind 
for mult i pole mag- 



netic field n > 5 ([Manchester fc Tavlorl 119771) . At 
the same time the measurements of pulsar frequency 
second derivatives v have shown that their values are 
much larger than expected for standard spin-down 
law and are even negative for about half of all pul- 
sars. The cor responding braking i ndice s range from 



10 b to 10" (D'A 



1993: Hobbs et al 



essandro et al. . 1993t ChukwudeL 
20041 ). 



It was found that the significant correlations 
between \v\ (\P\) and v (P) demonstrate the 
fact that the absolute values of the v and P 
are larger for younger (with large \v\) pulsars 



Preprint submitted to Elsevier 



1 February 2008 



(ICordes fc Downs! Il985t lArzoumanian et al . 1994 ; 
Lvnel . Il999h ." The anomalously high and negative 
values of v and n may be interpreted as a result 
of low- frequency components of the "timing noise" 
- a complex change of pulsars' rotation al phas e 



within a timespan (ID'Alessandro et all 119931 ) . 
Or, as a result of an y long-term influence o n the 



1977 



pulsars' spin-down. (IGullahorn fc Rankinl . 
Demiahsky &: Proszynskilll979h 

It is clear that the timespan of observations is 
by no means intrinsic to pulsar physics. Indeed, the 
variations of rotational parameters may take place 
on larger timescales as well. However, the timescale 
of observations naturally divides these variations 
into two separate classes of manifestations: (i) the 
well-known "timing noise" - the residuals in respect 
to the best fit for the timing solution and (ii) the 
"long term timing noise" - the systematic shift of 
the best fit coefficients (i.e. in the measured values 
of v, i>, v) relative to some mean or expected value 
from the model. 

Up to date we know nearly 200 pulsars for which 
the timespan of observations is greater than 20 
years, and the valu es of their i> still tu rn out to be 
anomalously large ( Hobbs et al . 2004f ). 



For example, for the PSR B1706-16 pulsar, varia- 



tions of v with an amplitude of 10" 



-24 



have been 



de tected on a timesc ale of several years (see Fig. 7 



Hobbs et al. (2004)), with the value of v depend- 



ing on the time interval selected. However, the fit 
over the entire 25 year timespan gives a value of v = 
3.8 • 10~ 25 s~ 3 with a few percent accuracy (which 
leads to a braking index w 2.7 • 10 3 ). 

In the current work we provide observational ev- 
idence of the non-monotonous evolution of pulsars 
on timescales larger than the typical contemporary 
timespan of observations (tens of years), using the 
statistical analysis of the measured v, v and v. We es- 
timate the main parameters of such long-term vari- 
ations and discuss their possible relation to the low- 
frequency terms of timing noise. We have also de- 
rived the parameters of pulsar secular spin-down. 
We plotted the v — v diagrams for 295 and for 1337 
"ordinary" radiopulsars. The smaller subset consists 
of pulsars with measured v. The bigger one - of all 
the "ordinary" pulsars wit h measured v and v. take n 



from the ATNF database (jManchester et alll2005T l 



We have found a good correlation beetween these 
two parameters and determnined the mean slope of 
the v — v distribution, which is in agreement with 
n ~ 5. 
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Fig. 1. The v — v diagra m for 295 pul s ars. T he figure shows 
the pulsars taken from ijfiobbs et all l2004h as circles, and 
the objects measured by other groups as squares. Open sym- 
bols represent the relatively young pulsars associated with 
supernova remnants. Analytical fits for both positive and 
negative branches are shown as solid lines. Measurment er- 
rors are shown as error bars. 



2. Statistical analysis of the ensemble of 
pulsars 

As was stated above, earlier works have shown the 
possibility of long-term variations of pulsars' rota- 
tional frequency. Our statistical analysis shows that 
the measurements of most the v reflect the pulsars' 
spin-down evolution on a timescale larger than the 
duration of observations, and uses the parameters of 
295 pul sars. From the 389 object s of the ATNF cat- 
alogue (jManchester et aL . 2005) with known v we 
compiled a list of "ordinary" radio pulsars with P > 
20 ms, P > 10 -17 s/s, and with relative accuracy 
of second derivative measurements better than 75%. 
We excluded recycled, anomalous and binary pul- 
sars. 26 supplement a ry pu l sars from o t her s ources 
(|D'Alessandro et all 1 19931 : IChukwudd . Il993h were 
added. The parameters of all pulsars were plotted 
on the v — v diagram (Fig. 1). 

The basic result of the statistical analysis of this 
data is a significant correlation between v and v, 
both for 168 objects with v > (correlation coef- 
ficient r w 0.90) and for 127 objects with v < 
(r w 0.85). Both groups follow nearly linear laws, 
however they are not exactly symmetric relative to 
v = 0. We divided both branches into 6 intervals 
of i>, computed the mean values and their standard 
deviations of i>± for each interval. We rejected the 
hypothesis of the branches symmetry with a 0.04 
significance level. Also, the absolute values of i) + 
are systematically larger than the corresponding V- 
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Fig. 2. The obvious ii — dependence for pulsars with 
the measured second derivative. The open symbols represent 
pulsars associated with SNRs. The solid line represents the 
best fit, the dotted lines — 1-cr range. The index of power 
— 1.16 ± 0.02 is significantly different from —1.0 which indi- 
cates the presence of the v — v correlation. 

(the difference is positive in 5 intervals out of 6) and 
the difference of analytical fits to branches is pos- 
itive over the — 10~ n -. — 10~ 15 s~ 2 interval of v. 
These are the arguments in favour of a small posi- 
tive assymetry of the branches. 

We found an obvious correlation of v with the 
characteristic age t c h = (r = 0.96, see Fig. 2): 



£ _ _2Q-6.53±0.12 T -1.16±0.02 



(1) 



This correlation will be discussed below in the light 
of the v — v dependency. The main feature of (1) 
is that the index of power —1.16 ± 0.02 is signifi- 
cantly different from —1.0. It is obvious that v would 



be ~ T~ h if v and v were completely uncorrelatcd. 
Therefore, the derived result provides an argument 
in favour of the v — v correlation. 

The v and r c h are nearly proportional, which leads 
to a significant correlation of T c h both with i) (r = 
0.85 for the positive branch and r = 0.75 for the 
negative one, Fig. 3) and with n (r = 0.75 and r = 
0.76 correspondingly, Fig. 8). 

The correlations found are fully consistent with 
the results published in dCordes k. Downs! . 1985 : 
Arzoumanian et al. , 19941 : lLvnei . ll999l ). as well as in 
( Urama et all 2*006| ). However, the branches with 
v > and i) < in those works were not analysed 
separately from each other (not as \i/\). 

Young pulsars confidently associated with super- 
nova remnants are systematically shifted to the left 
in Fig. 1 (open symbols). The order of magnitude 
of their physical ages roughly corresponds to that 
of their characteristic ages. This means that any de- 
pendence on v or T c h reflects the dependence on pul- 
sar age. The v — i> diagram (Fig. 1) may be inter- 



Fig. 3. The v — Tch dependence for pulsars with the measured 
second derivative. The open symbols represent pulsars asso- 
ciated with SNRs. This diagram is an argument in favour 
of the idea that pulsars evolve from large to small \v\ values 
independently from the i> sign. 



preted as an evolutionary one. In other words, each 
pulsar during its evolution moves along the branches 
of this diagram while increasing the value of its v 
(which corresponds to the increase of its characteris- 
tic age). However, there is an obvious contradiction: 
on the negative branch, v, being negative, may only 
decrease with time (since v is formally the deriva- 
tive of v), and the motion along the negative branch 
may only be backward! This contradiction is easily 
solved by assuming non-monotonous behaviour of 
i'(t), which has an irregular component (6v) along 
with the monotonous one (p ev ) , where the subscript 
"ev" marks the evolutionary value. 

The characteristic timescale T of such variations 
must be much shorter than the pulsar life time and 
at the same time much larger than the timescale 
of the observations. As it evolves, a pulsar repeat- 
edly changes sign of i>, in a spiral-like motion from 
branch to branch, and spends roughly half its life- 
time on each one. The asymmetry of the branches 
reflects the positive sign of iy e v(t), and therefore, sec- 
ular increase of i> ev (t) (i.e. all pulsars in their sec- 
ular evolution move to the right on the v — v di- 
agram). The systematic decreasing of the branches 
separation reflects the decreasing of the variations 
amplitude and/or the increasing of its characteristic 
timescale. 

Any well known non-monotonous variations of 
v(t), like glitches, microglitches, timing noise or pre- 
cession, will manifest themselves in a similar way on 
the v — v diagram and lead to extremely high va l- 
ues of v ( Shemar fc Lvne . 19961 : IStairs et al.l . l2000h . 
However, their characteristic timescales vary from 



weeks to years, and they are detected immediately. 
But here the variations on much larger timcscalcs 
are discussed, and their study is possible only sta- 
tistically, assuming the ergodic behaviour of the 
ensemble of pulsars. 

In addition, the apparent branches separation on 
the i> — v diagram is only due to the logarithmic scale 
of the plot. The spread of v values in each branch 
reaches 3 orders of magnitude, i.e. the pulsars on 
the diagram cover almost fully the range of possible 
v values (for each value of v) . Moreover, we observe 
the lack of pulsars near i) ~ because (i) the present 
accuracy of the v measurements is no better than 
10 -29 s~ 3 and (ii) in this area pulsars move faster 
than anywhere, because here v w const, v changes 
its sign and hence the third derivative of v has an 
extremum. 

3. Non-monotonous variations of pulsar 
spin-down rate on large timescales 

Variations of the pulsar rotational frequency may 
be complicated - periodic, quasi-periodic, or com- 
pletely stochastic. Generally, it may be described as 
a superposition 

v(t) = u ev (t)+5v(t), (2) 

where v ev (t) describes the secular evolution of pul- 
sar parameters and 5v(t) corresponds to the irregu- 
lar variations. Similar expressions describe the evo- 
lution of v and v after a differentiation. 5u{t) satis- 
fies the obvious condition of zero mean value: 

< Si)(t) > t = (3) 

over the timespans larger than the characteristic 
timescale of the variations. The amplitude of the ob- 
served variations of v is related to the dispersion of 
this process as 

V6i> =A i> = yf< (Si)) 2 > (4) 

The second derivative values on the upper v+ and 
lower i>_ branches in Fig. 1 may be approximately 
described as 

i> ± (t) = i> ev (t)±A i> (t) (5) 

for each pulsar. This equation describes an "aver- 
age" pulsar, while the spread of points inside the 
branches reflects the variations of individual param- 
eters over the pulsar ensemble and reaches 4 orders 
of magnitude. 
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Fig. 4. The v — v diagram for the pulsars with the measured 
second derivative. The filled symbols are objects with pos- 
itive v, the open ones - with negative v. The behaviour of 
both subsets is the same. The solid line represents the best 
fit corresponding to the braking index of n ~ 5, the dotted 
lines represent the l-cr range. 

The second derivative v is the only parameter sig- 
nificantly influenced by the timing variations (see 
Section 4). Thus one can assume that the measured 
values of v and v may be considered to be evolution- 
ary, u ev and u ev (since 8v and bv are small). 

3.1. v — v diagram 

Using the relations described above, the secular 
behaviour v(t) (or, v(y)) may be found by plotting 
the studied pulsar group onto the v~v diagram (Fig. 
4). The objects with v > and i) < are marked as 
filled and open circles, correspondingly. It is easily 
seen that the behaviour of these two sub-groups is 
the same, which is in agreement with the smallness 
of the pulsar frequency variations in respect to the 
intrinsic scatter of v(y). However, a strong correla- 
tion between v and v (r w 0.7) is seen, and 

= -Cv n , (6) 

where C = 10" 15 - 26±1 - 38 and n = 5.15 ± 0.34. So, 
the secular evolution of an "average" pulsar goes ac- 
cording to the "standard" spin-down law with n ~ 
5! This result is very interesting on its own, espe- 
cially since the v and v are always measured as in- 
dependent values. 

Note that the width of the fit on Fig. 4 is quite 
small - only about 1.5 orders of magnitude. If the 
spin-down is even approximately close to being de- 
scribed by the vacuum dipole model, then C should 
scale as (B sin x) 2 , where B is the polar field and 
X is the magnetic inclination angle. But the range of 
(Bq sin%) 2 over the pulsar population is expected to 
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be of many orders of magnitude. Therefore, Fig. 4 
once again shows that a simple vacuum dipole model 
is not adequate to observations. 

It is important that there is no contradiction be- 
tween the derived value of n ~ 5 and values of the 
braking index derived in the usual manner using v 
(these values lie in the range from approximately 
-10 6 to 10 6 ). The measured v and v are weakly af- 
fected by the timing noise, so the pulsars move along 
the v — v distribution mostly according to the sec- 
ular component of v. At the same time the combi- 
nation vv /v 2 = n includes both the secular and ir- 
regular components of i), which leads to anomalous 
braking indices. 

The braking index n w 5 allows us to describe 
the pulsars spin-down as "quadrupole-like" . But it 
is unlikely that such spin-down is due to the sim- 
ple quadrupole structure of the NS magnetic field 
and the corresponding radiation. Indeed, in case of 
vacuum approximation, for the rotating quadrupole 
with magnetic field strength Bq(Rms) ~ 10 12 G 
(where Rms ~ 10 6 cm), the value of the param- 
eter C from Eq. (6) would be as small a s 10 
(jManchester fe Tavlorl . Il977t iKrolikl . Il99lf ). which 
is ten orders of magnitude smaller than the values 
measured for the 295 and 1337 pulsars (see Figures 
4 and 5) . Moreover, the total contribution of all the 
high order multipole components (quadrupole etc.) 
is much smaller than that of the dipole component. 

Howcwcr, it is clear that the described sim- 
ple quadrupole spin-down does not satisfy the 
modern concept of th e pulsars' magneto s pheres 
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not being vacuum dGoldreich fe Julian . 



Manchester fc Tavloilll977tlBeskin et al 
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1993). 



A number of papers were published in the last 
years, where the possibility of quadrupole compo- 
nents of pulsars ' magn etic fields was discussed (e.g. 
Zane fc Turollal (|2006h ). Thus, at the same time, 
there is little reason to fully reject the hypothesis of 
a more complex pulsars' magnetic fields structure 
than only dipolar. 

An additional argument in favour of the v — v cor- 
relation is, as has been stated above, the dependence 
between v and r c h, where the index of power is not 
equal to —1.0. Indeed, since t c h = — and if the 
v — v correlation is absent, then v may only bring 
in some scatter to the i> — r c h dependence without 
its slope changing (on a logarithmic scale). But, if 
v and v are really correlated, and the slope (brak- 
ing index) is about 5, the v — t c h dependence should 
show a slope of about —1.2 which is roughly consis- 
tent with the measured value (see Eq. 1). Precisely, 
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Fig. 5. The v — u diagram for the 1337 pulsars from the ATNF 
database. The distribution shows quite a strong correlation 
(r = 0.60). The slope of the distribution corresponds to 
n = 5.48 ± 0.18. The slope of the dash-dot line on the plot 
corresponds to n = 2.5. It is clear that if the values of n for 
young pulsars are close to 2.5-^3.0, then they should change 
to higher values as the pulsars evolve. 

the slope of the i> — r fit, a, and n are related by the 
equation: 



(7) 



a - 1 



So it clearly seen that n is strongly dependent on a. 
The value of n derived from a = 1.16 ±0.02 is 7.25 ± 
0.78. It differs from 5.15±0.34 on a 2.5<r significance 
level. This value is less than the standard 3a but 
close to it. 

In a framework of pulsar spin-down analysis 
we ploted the v — v diagram for the 1337 pulsars 
(Fig. 5) taken from the ATNF pulsar database 
( Manchester et ail 2005f h This subset does not 
include recycled, binary and anomalous pulsars, 
and satisfies the criteria described in section (2). 
The braking index derived from the plot is n = 
5.48 ± 0.18. This value is in a good agreement with 
the result described above. The precision of the 
slope measurement on Fig. 5 is so high mostly due 
to the large number of pulsars used. 

The mean slope of the fit for the v — v distribution 
represents the average value of n for pulsars during 
their lifetime (« 5). Up to date there are several 
young pulsars with measured braking indices that 
are both accurate and precise. These are all smaller 
than but close to 3. The fact that the measured value 
of n was found to be greater than 3, could mean that 
pulsars do not evolve with a constant n (< 3). 

Indeed, if the pulsars are really born with n < 3, 
and assuming they do not change this value dur- 
ing their evolution, then as they grow older, they 
should appear below the main distribution on the 
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v diagram: approximately in the area with v 



10 



-16 



10 



-15 



s 2 and v ~ 0.1 Hz. 



As an example, the line with a slope which cor- 
responds to n = 2.5 is shown on Fig. 5. This line 
significantly deviates from the main distribution in 
the area of the oldest pulsars. If the pulsars were to 
evolve along this line, they would arrive in the area 
below the main distribution as they aged. 

Thus, the pulsars start their evolution with n < 3, 
which seems to be typical, and are likely to change 
their n to higher values as they evolve. Therefore, 
the authors consider to be reasonable the idea of 
young and old pulsars evolving with different values 
of n. These values are smaller for younger pulsars 
and larger for older ones. 

Such changing of n can be explained, for exam- 
ple, by the changing of the dominating spin-down 
mechanism. In general, a number of scenarios were 
proposed, where pulsars' braki ng index would in- 
crease with time (see for example lCordes fc ChernofI 
1998 ) and references therein; see also iRuderman 



2006)) 



On the other hand, the increasing of the braking 
index may be explained by the evolution of the pa- 
rameter C from Eq. (6). In case of a simple dipole 
spin-down C cx (Bq sin%) 2 and it is not constant for 
pulsars during their lifetimes. The magnetic field de- 
cay a nd magnetic inclination a ngle evolution (x < 



0, see Davis fc Goldstein (Il970h ) will lead to the dc 



creasing of C for the old pulsars. Hence, the abso- 
lute values of v will be systematically smaller than 
the "unshifted" ones (when C = const). Thus the 
slope of the v — i> distribution (braking index) will 
increase. 

A similar situation also tak es place in the mode l 
of electric current spin-down (jBeskin et all 1993). 
In this case, the relation v = —Kv 3 is also true, but 
K cx (Bo cos x) 2 with \ > 0, instead of (Bq sin x) 2 
with x < 0. 

In general, any pulsars' spin-down law should de- 
pend on the polar magnetic field Bq and, very likely, 
on the magnetic inclination angle x- Therefore, the 
idea described above should remain valid. If so, then 
the measured n rj 5 means that pulsars' intrinsic 
braking index may be less than 5. 

The provided analysis strongly suggests the pres- 
ence of a v — v correlation for ordinary pulsars with 
a mean slope (braking index) close to 5. Younger 
pulsars seem to evolve with a lower n than the older 
ones. The pulsars' magnetic field may have a more 
complex structure than dipolar, however there is lit- 
tle reason to insist that n ~ 5 is due to a simple 
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Fig. 6. The v — v diagram with the simple variations model. 
The solid line is the amplitude Ao = a v of the frequency 
second derivative variations, the dashed line is the secular 
term i>e V , and the dot-dashed lines are the envelopes of the 
oscillations v ev ± Oy with 1-cr ranges (dotted lines). Each 
pulsar spends the majority of its lifetime at or very near the 
envelopes. The pulsars are labeled as on Fig. 1 

quadrupole spin-down. It may be a result of the pul- 
sars' magnetic field decay and/or the evolution of 
the magnetic inclination angle. 

3.2. Some numerical results 

From Eq. (6) we may easily determine the relation 
between i^„„ and i> as 



v ev = nC™ (-z>) 2 



(8) 



which is shown in Fig. 6 as a thick dashed line. The 
same relation may be also estimated directly by us- 
ing the asymmetry of the branches seen on Fig. 1 as 

1 



(9) 



where i)± are defined in Fig. 1. Such estimation, 
while being very noisy, is positive in the — 10 -11 ~ 
— 10~ 15 s~ 2 range and agrees quantitatively with the 
previous one. 

The amplitude of the v oscillations, may be 
easily computed in a similar way, by using i) + and 
i>_, as = \{i> + - £>_). 

The behaviour of pulsars according to the derived 
relations is shown in Fig. 6. This simple variations 
model describes the observed branches, both posi- 
tive and negative, rather well. We interpret the ab- 
sence of negative branch objects with i> < —10 
s~ 2 , i.e. with r c h < 10 4 , as a prevalence of the sec- 
ond derivative's secular component over the varying 
one (An < v ev ) in this region. Older pulsars begin to 
change the sign of v because of spin rate variations. 
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Fig. 7. The v — v diagram for low noise pulsars. This di- 
agram shows a possible difference between long timcscale 
variations, discussed in the current work and the well-known 
"timing-noise" , because some of the pulsars (especially with 
low \i>\) still have anomalous braking indices. 

4. Discussion and conclusions 

In general, it is impossible to estimate the ampli- 
tudes of the frequency and its first derivative varia- 
tions A v and Ay from the amplitude of the second 
derivative only (the knowledge of its complete power 
density spectrum is needed). However, if the spec- 
tral density is relatively localized and some charac- 
teristic timescale T of the variations exists, it is pos- 
sible to set some limits on it. A rough estimation 
is A v ~ Ai,T, Ac ~ A V T and A v ~ A V T 2 . On a 
long timescale the variations can not lead to pulsar 
spin-up, and therefore, the variations of frequency 
first derivatives are much smaller than their secu- 
lar values, and Ay ~ A V T <C 0, so T <C v/A v . So, 
for the -10" 12 < v < -10" 
sponding values of A v from 10 _2,; 
the characteristic timescale T„ 



15 s 2 range and Carre- 
s' 3 to 10" 26 s" 3 , 
-up - 10 11 s. Also, this 
characteristic timescale is obviously larger than the 
timespan of observations: 50 < T < 3 • 10 3 years. 
Assuming the constancy of T during the pulsar evo- 
lution and therefore the change of A v with time, we 
get A v ~ 10 -3 -7- 10~ 7 Hz. For such a model the pul- 
sar frequency varies with the characteristic time of 
several hundred years and the amplitude from 10~ 3 
Hz for young objects to 10~ 7 Hz for older ones. 

The characteristic timescale T measured by the 
method above does not depend on the derived value 
of n ~ 5. Therefore, even if the braking index is 
significantly bigger than 5, the values of T and A v 
will not be changed. 

The physical reasons of the discussed non- 
monotonous variations of the pulsar spin-down rate 



Fig. 8. The n — T jj diagram. The older pulsars have a larger 
spread of braking index values. But some of them show a 
very low timing noise within the timespan of observation 
(see Fig. 7). 

may be similar to the ones of the timing noise on 
a short timescale. Several pr ocesses had been pro- 
posed for their explanation (jCordes fc Greenstein . 
19811 ) - from the collective effects in the neutron 
star superfluid core to the electric current fluctua- 
tions in the pulsar magnetosphere. Whether these 
processes are able to produce long timescale varia- 
tions is yet to be analysed. On a short timescale, the 
pulsars show different timing behaviour. But on the 
long timescale their behaviour seems to be alike. 

The argument in favour of the similarity between 
the discussed variations and the timing noise is the 
coincidence of the timing noise v amplitude ex- 
trapolated ac c ordin g to its power spectrum slope 
(jBavkal et all I1999T ) to the timescale of hundreds 
of years, with the Ay derived from our analysis for 
the same v, i.e. the same ages. 

At the same time, there are several low noise 
pulsars with large or negati ve v (see Fig s . 7 an d 8). 
For 55 pulsars studied in (|Hobbs et all , 12004) the 
timing noise is nearly absent (RMS < 1 • 10 _3 P). 
Some of them have anomalous values of v , which 
are well consistent with the \v\ - v correlation 
( Cordes fc Downsl 119851 lArzoumanian et all . Il994h 
and have a wide range of v. This shows a possi- 
ble difference between timing noise and the long 
timescale variations described above. 

In any case, the principal point is that all the 
pulsars evolve with long-term variations, and the 
timescale of such variations significantly exceeds 
several tens of years. That explains the anomalous 
values of the observed v and braking indices and 
gives reasonable values of the underlying secular 
spin-down parameters. 
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